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Abstract. We present the results of a Near-Infrared deep photometric survey of a sample of six embedded star clusters in 
the Vela-D molecular cloud, all associated with luminous (~ 10^ Lq) IRAS sources. The clusters are unlikely to be older 
than a few 10* yrs, since all are still associated with molecular gas. We employed the fact that all clusters lie at the same 
distance and were observed with the same instrumental setting to derive their properties in a consistent way, being affected by 
the same instrumental and observational biases. We extracted the clusters' K Luminosity Functions and developed a simple 
method to correct them for extinction, based on colour-magnitude diagrams. The reliability of the method has been tested by 
constructing synthetic clusters from theoretical tracks for pre-main sequence stars and a standard Initial Mass Function. The 
clusters' Initial Mass Functions have been derived from the dereddened K Luminosity Functions by adopting a set of pre-main 
sequence evolutionary tracks and assuming coeval star formation. All clusters are small (~ 100 members) and compact (radius 
~ 0. 1 - 0.2 pc); their most massive stars are intermediate-mass (~ 2 - 10 Mq) ones. The dereddened K Luminosity Functions 
are likely to arise from the same distribution, suggesting that the selected clusters have quite similar Initial Mass Functions and 
star formation histories. The Initial Mass Functions are consistent with those derived for field stars and clusters. Adding them 
together we found that the "global" Initial Mass Function appears steeper at the high-mass end and exhibits a drop-off at ~ 10 
M0. In fact, a standard Initial Mass Function would predict a star with M > 22.5 Mq within one of the clusters, which is not 
found. Hence, either high-mass stars need larger clusters to be formed, or the Initial Mass Function of the single clusters is 
steeper at the high-mass end because of the physical conditions in the parental gas. 

Key words. Stars: formation - Stars: pre-main sequence ~ ISM: individual objects: Vela Molecular Ridge - Infrared: stars 



1. Introduction 

It is now well established by means of direct and indirect 
observations that most, if not all, stars are formed in groups 
rather than in isolation (Clarke et al. 2000 Lada & Lada 2003). 
Adams & Myers (2001) suggested that most stars form in small 
clusters with ~ 10-100 members. However, by studying a large 
sample of embedded clusters within 2 kpc, Lada & Lada (2003) 
concluded that 90 % of stars that form in embedded clusters oc- 
cur in rich clusters of > 100 members, but only < 4 - 7 % of 
embedded clusters emerge from molecular clouds to become 
bound clusters of Pleiades age. 

An important result that strongly constrains theories of 
massive stars and stellar clusters formation is that the stellar 
density of young stellar clusters seems to depend on the most 
massive star in the cluster. In low-mass star forming regions, 
stars are usually found to form in loose groups with typical den- 
sities of a few stars per cubic parsec (Gomez et al. 1993 1, while 
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high-mass stars are found in dense clusters of up to 10'* stars per 
cubic parsec, along with a large fraction of low-mass stars (e.g. 
the Orion Nebula Cluster, Hillenbrand & Hartmann .1998) . The 
transition between these two modes of formation should occur 
in the intermediate-mass regime, namely 2^ M/M© $15. 

In order to probe this transition, Testi et al. <1999> car- 
ried out an extensive Near-Infrared (NIR) survey searching for 
young clusters around optically visible intermediate-mass stars 
(Herbig Ae/Be stars) in the northern hemisphere. The main re- 
sult of this survey is that there is a strong correlation between 
the spectral type of the Herbig Ae/Be stars and the membership 
number of the stellar groups around them. Furthermore, there 
is compelling evidence that the most massive stars in their sam- 
ple are surrounded by denser, not only more populous, clusters. 
These findings are in qualitative agreement with models that 
suggest a causal relationship between the birth of a massive 
star and the presence of rich stellar clusters. The observed cor- 
relation and scatter, however, could also be explained in terms 
of random assembling clusters with membership size distribu- 
tion of the form g(N) ~ N^^-^ picking stars from a standard 
Initial Mass Function (IMF; BonneU & Clarke il999< . 
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As discussed in Testi et al. i2001 ^') there are two obser- 
vational strategies that may provide additional constraints on 
which of the two scenarios is the most appropriate: expand the 
sample of optically revealed young O and B stars to increase 
the statistic, and search for clusters in complete samples of lu- 
minous embedded sources in giant molecular clouds. Recently, 
de Wit et al. (.20041 12005> . following the first approach, in- 
vestigated the possible origin in clusters of field O-stars; they 
demonstrated that the majority of the so called "field" O-stars 
may have originated in clusters, but could not exclude that a 
small fraction of these are actually born in isolation. An exam- 
ple of the second approach is given in Massi et al. (2000, 2003) 
as described in Sect. |2] 

One fundamental property of clusters is their IMF; differ- 
ent methods to extract the IMF from photometric and/or spec- 
troscopic data are summarized in Lada & Lada (2003). These 
authors concluded that embedded young clusters exhibit strong 
similarities in their IMF with open clusters and field stars, al- 
though this "universal" IMF does not seem to characterize all 
star forming regions (Lada & Lada 2003). Scalo (1998), re- 
viewing all determinations of the IMF in different environ- 
ments, noted that either the observational uncertainties are so 
large as not to allow one to draw definite conclusions on the av- 
erage IMF and its variation, or the available data show strong 
evidence for IMF variations. This is a critical issue, since the 
IMF should represent the imprint on a stellar population of the 
physical conditions in the parental gas. 

Most of the main topics concerning the formation of star 
clusters are then still debated and ask for large observational 
efforts. In this paper, we report the results of a NIR study of 
a sample of small clusters in the giant molecular cloud "D" of 
the Vela Molecular Ridge. Small young embedded clusters can 
prove very interesting targets to address many of the illustrated 
points, since they 1) are associated with the intermediate-mass 
transition regime investigated by Testi et al. ([l999i), 2) are on 
the border between the different sizes envisaged by Adams & 
Myers (2001) and Lada & Lada (2003) as major star contribu- 
tors, and 3) do not host stars massive enough to heavily affect 
the parental environment, maintaining a stronger link with it. 
All imaged clusters are young and are associated with the same 
molecular cloud. On the one hand, this allows one to circum- 
vent many effects due to observational biases in comparing the 
results. On the other hand, local differences in the global IMF 
emerging from the molecular cloud can be searched for and in- 
vestigated. Furthermore, our data considerably improve on pre- 
vious NIR photometry of the same regions, both in sensitivity 
and in field of view. The Paper layout is as follows: in Sect.|2l 
the main features of the star forming region and the selected 
sample are summarized. Observations and data reduction are 
described in Sect. |3l The main results are reported in Sect. |3 
and discussed in Sect|5l Our conclusions are listed in Sect. |6l 

2. The Vela-D Luminous IRAS Sources 

The Vela Molecular Ridge (VMR) is a giant molecular cloud 
complex located in the outer Galaxy (257° < I < 274°) within 
the galactic plane (-5° < b < 5°). It was first mapped in the 
CO(l-O) transition (with low resolution) by Murphy & May 



Table 1. Selected fields. The last Col. lists the designation 
adopted by Liseau et al. (1991). 
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(1991), who divided the emission area into 4 main regions 
(named A, B, C and D) corresponding to local maxima. The 
kinematical distance is 1-2 kpc. 

Liseau et al. (1991) and Lorenzetti et al. (1993) studied the 
star formation in the VMR by selecting the IRAS point sources 
with red colours and combining their Far-Infrared (FIR) fluxes 
with NIR and 1 .3-mm photometry. They list a sample of IRAS 
sources which, based on the SEDs from NIR to FIR wave- 
lengths and the radial velocity of the associated molecular gas, 
may be considered as intermediate-mass (M ~ 2 - 10 Mq) 
analogues of Class I sources belonging to the VMR. Liseau et 
al. (1991) also discuss the issue of distance, concluding that 
clouds A, C and D are at li = 700 + 200 pc, whereas cloud 
B is further out (~ 2 kpc). Massi et al. (2000) analyzed JHK 
images of the fields towards 12 sources of the sample, namely 
those associated with cloud D, whose bolometric luminosities 
range from ~ 10^ to ~ 6 x 10^ L© (Massi et al. 1999). These 
authors found young embedded star clusters towards the IRAS 
sources with Lboi ~ 10^ L©. Five of these embedded clusters 
have been selected for the present work. In addition, we have 
included IRAS 08438-4340 (IRS 16, as designated by Liseau 
et al. 1991), a FIR source associated with both an embedded 
star cluster (Massi et al. 2003) and an HII region, close (in pro- 
jection) to the other ones. The six fields are listed in Tabled 
along with their equatorial coordinates and the nomenclature 
used by Liseau et al. (1991). 

3. Observations and data reduction 

The six fields towards the selected IRAS sources in the Vela 
D molecular cloud were observed at the ESO-NTT telescope 
on February 18, March 5, March 20 and September 25, 2000. 
Near-Infrared images through the standard J, H, and broad 
band filters were acquired using the Sofl instrument, with a 
plate scale of ~0.282"/pixel and an instantaneous field of view 
of ~5'. For each target we collected a set of dithered exposures 
for a total integration time of 15 minutes in each filter. The raw 
images were crosstalk corrected, flat fielded, sky subtracted, 
aligned and mosaiced using the special procedures developed 
for Sofl and standard routines within the IRAF package. 

The final image quality of our images was ~ 0.85" for all 
but the field surrounding IRS 16 which was observed under 
slightly worst seeing conditions (1.1"). Photometric calibra- 
tion was ensured by observations of a number of standard stars 
from the hst of Persson et al. (.1998^ . Source detection and aper- 
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ture photometry were performed using standard routines within 
IRAF. 

We compared our photometry with that of Massi et al. 
(1999) and Massi et al. (2003) to check for consistency. The 
latter was obtained from lrac2 JHK images. The mean differ- 
ences (in the sense lrac2 minus Sofl) are: ~ to -0.07 mag 
(excepted IRS 19 and IRS 21, with -0.34 and -0.16 mag) in 
K, -0.06 to 0.01 mag (excepted IRS 17 and IRS 21, with 0.3 
and -0.2 mag) in H and -0.1 to 0.08 mag (excepted IRS 17 
and IRS 16, with 0.26 and -0.16 mag) in J. Part of the largest 
mean differences may be attributed to problems with the aper- 
ture correction of some Irac2 data sets, as already noted by 
those authors (Massi et al. 2000; see also their Fig. 1 for IRS 17 
and IRS 19). The scatter of the differences around the average 
values are 0.11 to 0.26 mag in K, 0.12 to 0.19 mag (excepted 
IRS 16 and IRS 17, with 0.34 and 0.33 mag) in H and 0.21 
to 0.25 mag (excepted IRS 16 and IRS 17, with 0.36 and 0.30 
mag) in J. We found also a colour term in K of the order of 

0.2 X (H-K), certainly due to the different K (Irac2) and 

(Sofl) filters, that may account for another part of average and 
rm.s. differences. In H and J the colour term, if any, is much 
smaller: ~ 0.04 - 0.07 x (7 - H). We believe that the relatively 
large scatter in the magnitude differences may be explained by 
the differences in PSF sampling and sky estimate between the 
data sets and the intrinsic variability of young stars. Anyway, 
both data sets appear to be mutually consistent. 

4. Results 

In Fig.^we show the NTT/SofI A", images of the six surveyed 
fields. True-colour images (combining the images in all three 
bands) are shown in Testi et al. (2001b; see their Fig. 1). From 
those images it is immediately clear that we detected groups of 
very red sources in every fields. 

4.1. Mass sensitivity 

We estimate a completeness magnitude /T^ ~ 18 for all fields, 
excepted IRS 21 for which K'^^ ~ 17 (see Fig. 2 in Testi et al. 
2001b). Assuming a distance modulus of 9.22 mag, this allows 
an estimate of the minimum mass, as a function of age and red- 
dening, for a star whose probability to be found is still almost 
1 throughout each frame. Using the evolutionary tracks for 
pms stars given by Palla & Stabler (1999), the mass complete- 
ness limit is Mcompi ~ 0.1 Mq for Ay = 30 mag, /T, ~ 17 and 
an age of 10^ yrs and increases to ~ 0.4 M0 for an age of 10^ 
yrs. At this age, for Ay - 20 mag or = 18, Mcompi is still as 
low as 0.2 Mq. We find the same using the evolutionary tracks 
from D'Antona & MazzitelH (1994) for lO*" yrs old stars, but 
Mcompi ranges from ~ 0.5 Mq (Ay - 30 mag) to ~ 0.3 M0 
(Ay - 20 mag or K'^^ - 18) for 10^ yrs old stars. From Baraffe 
et al. (1998) the estimated Mcompi is ~ 0.15 Mq for 2 x 10*= yrs 
old stars and ^ 0.6 M0 for 10^ yrs old stars. It is clear that our 
images probe the clusters stellar population down to the hydro- 
gen burning limit for ages < 10^ yrs and as low as 0.5 Mq for 
an older population (10^ yrs) ; they are probably deep enough to 
probe the brown dwarf regime for ages ~ lO'' yrs. However, we 
caution against the decrease in completeness magnitude that is 



Table 2. Clusters' properties. The diameter is listed in Col. 2, 
the peak surface stellar density in Col. 3, the richness indicator 
(see text) in Col. 4 and the volume stellar density in Col. 5. 
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likely to occur towards the clusters' centres because of crowd- 
edness, diffuse emission and, sometimes, the PSF wings of the 
brightest stars. 

4.2. Clustering 

Although the star clusters are clearly visible in the images of 
Fig.^ we employ the same graphical representation as used by 
Massi et al. (2000) in order to better evidence their presence. 
Hence, we counted all Ks sources up to /T^ - 18 in ~ 30" x 30" 
squares ("bins"), at offsets of ~ 15" both in right ascension 
and in declination. Following Massi et al. (2000), we estimated 
the field star density (the "sky") and its 1 cr fluctuation by gath- 
ering all the bins from the 6 fields and constructing a global 
histogram of the number of bins with given counts. The ris- 
ing part of the histogram is well fitted by a Poissonian curve 
with a mean of 7.5 counts per bin (32 stars arcmin"^), so we 
assume that this represents the statistic of field stars. The ac- 
tual data distribution exceeds the Poissonian curve in the wing 
and, as noted by Massi et al. (2000), this indicates that the clus- 
tering is not due to sky fluctuations. The Poissonian fit yields 
a o" ~ 12 stars arcmin"^. Since our imaged fields are large 
enough, we constructed a histogram for each frame and by the 
same technique we obtained sky means ranging from 26 to 34 
stars arcmin"^ (and cr ~ 10-12 stars arcmin"^). These values 
compare quite well with those found by Massi et al. (2000), i. 
e., a sky mean ~ 20 stars arcmin"^ and a cr ~ 12 stars arcmin"^, 
considering that their photometry has a limiting magnitude 
~ 18 and, hence, a lower completeness magnitude than ours. 
Figure 121 shows the contour maps of surface star density from 
the Ks counting; contours start at 50 stars arcmin"^ (roughly 
the sky mean plus 2cr) in steps of 20 stars arcmin"^ (~ 2 cr). 
Clearly, all clusters have maximum densities well above the 
sky mean plus 3cr. This confirms that the structures found by 
Massi et al. (2000) on smaller field-size images are actually 
young star clusters. The clusters are compact enough to fit the 
2' X 2' f.o.v. of the images used by those authors. However, the 
larger f.o.v. allowed by Sofl and our much deeper imaging ev- 
idence some substructuring, such as the small source grouping 
(which we name "B") north-east of the cluster found by Massi 
et al. towards IRS 17 (which we name "A"). 

We determined the radial distribution of surface star density 
as described in Testi et al. (1999), shown in Fig.|3] Counting 
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Fig. 2. Star surface density (arcmin ^) from the Ks images of the 6 fields. Contours are in steps of ~ 2cr (20 arcmin ^) from ~ 2cr 
above the mean surface density of field stars. Towards IRS 17, two nearby subclusters (A and B) are clearly present. 
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Fig. 3. Radial distribution of star surface density (arcmin ) from the Kg images of the 6 fields. 



is performed on the sources with at least a detection in Ks 
and up to - 17, which is the smaller completeness mag- 
nitude we obtained throughout the fields. A significant density 
increase towards r = pc is noticeable in all the six fields. 
This also allows an estimate of the surface density due to field 
stars by averaging the radial density in the wing of the distribu- 
tion. Obtained values range between 14 and 25 stars/arcmin^ 
and are consistent with the one of 32 stars/arcmin^ we in- 
fer above, considering the latter is obtained for sources up to 
Ks = 18. The richness indicator Ic (Testi et al. 1999) is de- 
rived by integrating the radial stellar surface density, once sub- 
tracted the background-foreground contribution as estimated 
from the wing of the distribution. This is listed in Table |2 
along with the statistical error due to the sky determination and 
the Poissonian fluctuations in the surface density of the cluster 
members. Fig ure|3l shows that all clusters have a rather similar 
size and allows to estimate a radius ~ 0.1 - 0.2 pc (~ 0.5 - 1 
arcmin) at the assumed distance (see Table|2j. 

4.3. Cluster star population 

The nature of the star population associated with the embedded 
clusters can be studied through colour-colour (J-H\s. H-K,) 
and colour-magnitude (K,, vs. H-K,) diagrams. Exploiting our 
field of view, we compare the colours of sources within 2 ar- 
cmin of the cluster centre with those of sources outside. This 



distance is slightly more than twice the clusters' radius; in fact, 
we expect the radial profiles to exhibit a lower density tail far 
from the centre, so we chose a value of 2 arcmin in order to 
account for this. Radial density profiles on a much larger scale 
than our images, obtained from 2MASS data, actually confirm 
the presence of a tail up to at least twice the clusters' radius. 
Hence, the sky area outside this distance is less contaminated 
by clusters' members and more representative of the popula- 
tion of field stars. Figure |4] shows the colour-colour diagrams 
(CCDs) for objects within the chosen radius, whereas Fig.|5] 
shows the same but for objects outside, both up to K, = 18. 
The locus of main sequence (from Koornneef 1983) is drawn 
as a solid line and the dashed lines are parallel to the reddening 
vector (with crosses at intervals of Ay - 10 mag) according to 
the extinction law given by Rieke & Lebofsky (1985). The NIR 
sources > 2 arcmin far from the clusters' centres are essentially 
reddened stars, up to Ay ~ 10 - 15 mag; the dispersion around 
the main sequence and its reddening band denotes a large frac- 
tion of faint objects that may represent background field stars. 
There are few NIR sources with an intrinsic colour excess. 
Conversely, within the chosen radius, a different stellar pop- 
ulation emerges from the previous one; the objects now spread 
over a larger extinction range, up to Ay ~ 20 - 30 mag, there 
is less dispersion around the main sequence and the fraction 
of NIR sources with an intrinsic colour excess is much larger. 
This indicates that the peak extinction is greater towards the 
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Fig. 4. Colour-colour diagrams (i. e., J - H vs. H - Ks) of the 6 fields, for objects lying within 2 arcmin (i. e., roughly twice the 
clusters' radius) from the clusters' centres and up to Kg = 18. The solid line marks the main sequence (Koornneef 1983) whereas 
the dashed ones indicate the reddening law with crosses at intervals of Ay = 10 mag (Rieke & Lebofsky 1985). For the sake of 
clarity, an arrow corresponding to an extinction Ay = 10 mag has been also drawn, in the upper left panel. 



clusters' centres (hence, most background stars are reddened 
out) and that a population of very young stars is superposed to 
that of field stars which characterizes the outer edge of the im- 
aged fields. Including only sources with < 15 further limits 
the dispersion around the main sequence and reddening bands, 
but does not alter the spread along the extinction vector and the 
high fraction of NIR objects with intrinsic colour excess near 
to the clusters' centre. Also, note that the observed reddening 
follows quite well the drawn loci, confirming that the extinc- 
tion law towards the 6 fields is consistent with the standard one 
as given by Rieke & Lebofsky (1985). 

All the previous findings can be checked using colour- 
magnitude diagrams (CMDs). These are presented in Fig. |5] 
and Fig.0for all fields. The solid line is the locus of Zero Age 
Main Sequence (ZAMS) at the distance of Vela (700 pc) using 
the absolute magnitudes given by Allen (1976) and the colours 
given by Koornneef (1983). The location of BO, AO, KO and M5 
stars are labelled. The arrow is for a reddening Ay = 20 mag 
according to Rieke & Lebofsky (1985). It can be seen that the 
sources far from the clusters' centres (see Fig.Q are on average 
closer to the ZAMS (i. e., less extincted) and a large fraction of 
them is fainter than ~ 16 (i. e., they are mostly background 
field stars as noted above). Actually, the fraction of sources (up 
to Ks - 18) with Ks < 16 within/outside 2 arcmin from the 
centres ranges from 0.46/0.33 (IRS 16) to 0.36/0.29 (IRS 20). 



Furthermore, the sources within 2 arcmin from the clusters' 
centres are in a larger fraction far from the ZAMS, because of 
both their greater extinction and their intrinsic colour excess: 
the fraction of sources (up to Ks = 18) with H - Ks > 1.5 mag 
within/outside 2 arcmin from the centres ranges from 0.18/0.03 
(IRS 16) to 0.07/0.007 (IRS 18). This confirms that the clus- 
ters' stellar population is represented by very young, embed- 
ded stars, mostly located within 2 arcmin (or even less) from 
the centres and that it is indeed diff'erent from the one of field 
stars. 

4.4. The clusters' K Luminosity Functions 

The K Luminosity Function (KLF) of an embedded cluster is 
a simple tool to probe its Initial Mass Function (IMF) and age. 
As discussed, e. g., by Zinnecker et al. (1993), the KLF de- 
pends not only on the underlying IMF, but also on the mass- 
luminosity relation. For pms stars the mass-luminosity relation 
is a function of age. Hence, given the IMF, a cluster's KLF 
evolves with time and also changes according to the star for- 
mation history of the cluster. However, Muench et al. (2000) 
show that the KLF of a young cluster is much more sensitive to 
variations in its underlying IMF than to the star formation his- 
tory or the mass-luminosity relation. Also, a cluster KLF shifts 
to fainter magnitudes and widen systematically with age. This 
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Fig. 5. Same as Fig|3] but for objects further than 2 arcmin from the clusters' centres. 



should allow to put some constraints on the age of a young 
clusters, as well. 

The KLFs obtained from our photometry are shown in 
Fig. 2 of Testi et al. (2001b). For each field, two KLFs are 
determined: one using sources close to the cluster centre and 
one using sources at the image edge. Clearly, the two sets of 
distributions differ in shape, confirming that the clusters host 
a star population different from the field one, as inferred in 
Sect, lis] on the basis of CCDs and CMDs. In particular, the 
outer KLFs increase steadily and peak at higher magnitudes 
than the "central" KLFs. This maximum reflects the complete- 
ness limit at Ks- Since at the edges of our images the num- 
ber of cluster members is small with respect to field stars, 
we chose to construct "control" KLFs using the photometry 
obtained in the outer areas themselves of our frames. These 
should better characterize the contribution of field stars to the 
raw KLFs than those obtained using relatively distant control 
regions. We adopted a radius of 2' from the cluster centre to 
define the border between "centre" and "edge". The cluster's 
centre is assumed to coincide with the NIR counterparts of 
the IRAS sources identified by Massi et al. (1999). The reason 
for our choice of the radius is explained in Sect. 14.31 actually 
a larger radius would assure a much less contamination, but 
would also decrease the outer area available for estimating the 
control KLF 

Thus, we constructed the two KLFs per field by counting 
the Ks sources in bins of 0.5 mag, within and outside the se- 
lected radius. After normalizing the outer KLFs to the same 



area of the central ones, the clusters' KLFs were finally ob- 
tained by subtracting the former from the latter in each field. 
Yet, in order to derive the clusters' IMFs from their KLFs, 
source extinction and NIR excess have to be accounted for. 
Extinction causes a twofold problem. First, even with so nearby 
control fields the sources at the centre and those at the edge of 
each image are reddened in different ways, as evidenced by the 
CCDs. Furthermore, the single clusters' members are differen- 
tially reddened and, because the span in Ay may be more than 
10 mag (i. e., more than 1 mag in Ak), this cannot be neglected 
even binning the KLFs at a large interval. Thus, we tested two 
methods in order to account for the effect of extinction on the 
KLFs: modelling the reddening of theoretical KLFs, and trying 
to "deredden" the observed KLFs. We found the latter method 
more satisfactory and we discuss it in the following. 

4.5. Modelling K-Luminosity functions 

We performed several tests with synthetic KLFs, which are 
briefly discussed in the following sections. To do that, we de- 
veloped a Montecarlo procedure to construct theoretical KLFs 
by assuming an underlying IMF, an age and a star formation 
history. The adopted mass-luminosity relation has been derived 
from the evolutionary tracks of Palla & Stabler (1999). The pro- 
cedure makes an estimate of the number of cluster members on 
the basis of the bolometric luminosities that is selected, using 
an empirical relation between Ic and Lboi derived from the data 
of Testi et al. (1998, 1999). The NIR colour excess has been 
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Fig. 6. Colour-magnitude diagrams (i. e., Ks vs. H-Ks) of the 6 fields, for objects lying within 2 arcmin from the clusters' centres 
and up to Ks - 18. The solid line is the Zero Age Main Sequence (from Allen 1976 and Koornneef 1983) at the distance of the 
VMR, whereas the arrow is for a reddening Ay - 20 mag (Rieke & Lebofsky 1985). The location of BO, AO,MO and M5 ZAMS 
stars is labelled. 



modelled as in Hillenbrand & Carpenter (2000) and it is also 
possible to add to the NIR fluxes a random reddening (accord- 
ing to a Gaussian statistic) superposed to a constant value. 

4.6. Dereddening the observed K Luminosity 
Functions 

There are indications that the surveyed clusters are at most 
a few 10^ yrs old: their compactness, the high spread in ex- 
tinction displayed by their members and a significant fraction 
of sources exhibiting a NIR infrared excess in the CCDs (see 
Fig-E}- According to Lada & Lada (2003), the embedded phase 
of cluster evolution lasts 2-3 Myrs and clusters older than 
5 Myrs are rarely associated with molecular gas. All our ob- 
served clusters are associated with molecular gas (see Massi et 
al. 2000, 2003, 2005). 

Their youth was also confirmed through a preliminary mod- 
elization, constructing reddened theoretical KLFs based on a 
rough estimates of the extinction distribution in the clusters ac- 
cording to the CCDs. The comparison with the observed cluster 
KLFs indicates that the clusters' age roughly ranges between 
~ 10^ and ~ 10^ yr This suggested us a simple method in or- 
der to deredden the observed KLF, which is illustrated in Fig.|8] 
Figure Ell details a CMD (K vs. H - K) including flie ZAMS 
locus (solid line) and three isochrones from the evolutionary 
pms tracks of Palla & Stabler (1999; dashed lines), all for a 



distance of 700 pc. Clearly, pms stars from 10^ years old to the 
ZAMS span quite a limited range in colour (0.2-0.4 mag) and, 
at the same time, the slope of the reddening law is small. This 
means that the stars may be dereddened on a CMD simply by 
projecting them back along the reddening vector onto a mean 
locus (thick line in Fig. |8} chosen to be located between the 
isochrones. The error is small; e. g., for pms stars older than 10* 
yr, the error in K due to the procedure is ^ 0.2 mag, less than 
the adopted KLF binning interval (0.5 mag). Figure [S^ shows 
the worst case, one of a 0.1-Mq 10^-yrs-old pms star, whose 
obtained K is 0.28 mag less than its actual K after deredden- 
ing to the mean locus. Hence, only if a large fraction of cluster 
members were very young low-mass stars, the dereddened KLF 
could be significantly deformed. Even so. Fig. [SJ? shows that 
this result depends on the adopted evolutionary tracks: using 
those of D'Antona & Mazzitelli (1994) somewhat reduces the 
discrepancy. Finally, Fig.lSj; shows that an error on the cluster 
distance cannot cause a significant deformation of the dered- 
dened KLF with respect to the actual "unreddened" one. 

How well the dereddening procedure works is illustrated in 
Fig.m A model "unreddened" KLF is derived from one of the 
statistical realizations of a cluster with the IMF of Scalo (1998), 
truncated at 0.1 Mq, and using the evolutionary tracks of Palla 
& Stabler (1999). For each star the JHK magnitudes are scaled 
to a distance of 700 pc and reddened by adding a random value 
(ranging from to 20 mag) to an overall Ay - 20 mag (screen). 
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Fig. 7. Same as Fig|6l but for objects more than 2 arcmin from the clusters' centres. 



assuming the extinction law of Rieke & Lebofsky (1985). The 
synthetic sample is then projected onto our adopted mean locus 
in the CMD and a "dereddened" KLF is constructed. Cases a 
and d in Fig.Ocompare unreddened and dereddened KLFs for 
a coeval cluster of 10^ years old and for one with continuous 
star formation from 10^ to 10^ years old, respectively, when no 
NIR excess is superposed to the photospheric emission. It is ev- 
ident that the "model" and the "dereddened" KLFs are almost 
identical. 

Next, we tested how much the dereddened KLFs are af- 
fected by the NIR infrared excess of the clusters' members. 
Thus, we added to each star of the synthetic sample a random 
K contribution modelled as in Hillenbrand & Carpenter (2000). 
As shown in Fig.|5]for cases b and e (same as a and d but with 
a NIR excess), the high luminosity end of the "dereddened" 
distribution stretches with respect to the "unreddened" KLF. A 
slightly better matching is obtained by shifting the mean locus 
in H-K redward of 0.4 mag (roughly half the maximum H-K 
excess found by Hillenbrand & Carpenter 2000), as illustrated 
by cases c and f. 

The effect of changing the distance of the synthetic cluster 
or the evolutionary tracks is shown in Fig. The cluster of 
case a above (coeval star formation, 10^ yrs old) was again syn- 
thesized but using the tracks of D'Antona & Mazzitelli (1994), 
both without (case a) and with (case d) NIR excess; it is evident 
that the choice of the evolutionary tracks underlying the mod- 
elled cluster does not affect much the dereddened KLF which 
is obtained, whereas the effect of the NIR excess is the same as 
inFig.|9l 



The mean locus has been derived "theoretically" by assum- 
ing a distance of 700 pc, so we tested also what could happen if 
the actual clusters were at a different distance. Thus the cluster 
of case a in Fig.|9l(coeval star formation, 10^ yrs old, tracks of 
Palla & Stabler 1999) was put at different distances. In Fig. 1101 
the dereddened KLF is shown for a distance of 500 pc (b with- 
out NIR excess, e with NIR excess) and 900 pc (c without NIR 
excess, f with NIR excess). Clearly, the "dereddened" KLF al- 
ways reproduces quite well the "unreddened" KLF when no 
NIR excess is included. When a random NIR excess is added, 
the deformation of the "dereddened" KLF at the high luminos- 
ity end does appear similar for all cluster models illustrated in 
Figs.|9land[l0l 

We repeated this kind of visual test for a number of different 
synthesized clusters, not shown here, varying the number of 
members, the age and the star formation history. We found that 
when no NIR excess is included the "dereddened" KLF always 
matches quite well the "model" one for coeval clusters of 10^ 
yrs old, whereas behaves slightly worse for 10^ yrs old clusters. 
Anyway, the surveyed clusters are very unlikely to be so old. 
The "dereddened" KLF for extremely young clusters (^ 10^ yrs 
old) should appear narrower than the actual one, mimicking an 
even younger age; we did not perform experiments on such an 
unevolved star population, since probably the obtained KLF is 
much more sensitive to the adopted evolutionary model (see, e. 
g.,Fig.|H}. 

We carried out also another simple test in order to get an 
unbiased assessment of the goodness of the dereddening tech- 
nique. We constructed 6 synthetic coeval clusters of age rang- 
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Fig. 9. Comparison of model "unreddened" KLFs (open squares, dashed line), from one of the statistical realizations of a cluster 
at 700 pc as explained in the text, and "dereddened" KLFs (solid line), obtained adding a random extinction to each member 
of the cluster and then dereddening it according to our proposed method. Case a is for a coeval 10^ years old cluster without 
NIR excess, case b includes a NIR excess and in case c the mean locus for dereddening has been shifted by 0.4 mag redward in 
H — K m order to partially compensate for the NIR excess. Cases d,e and/same as a, b and c, but for a cluster modelled with a 
continuous star formation from 10^ to 10^ years. Bins with a number of counts N - Q have been set to log = -1 and a constant 
equal to 7,4.5 and 2.0 has been added to log for ad, be and cf, respectively. 
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Fig. 10. Comparison of model "unreddened" KLFs (open squares, dashed line), for the same cluster (of 10^ yrs) as case a of Fig.|9j 
and "dereddened" KLFs (solid line), obtained adding a random extinction to each member of the cluster and then dereddening 
it according to our proposed method. Cases a and d are for the evolutionary tracks of D' Antona & Mazzitelli (1994) without {a) 
and with id) a random NIR excess. In cases b (without NIR excess) and e (with NIR excess) the cluster has been put to a distance 
of 500 pc, whereas in cases c (without NIR excess) and /(with NIR excess) it has been put to a distance of 900 pc. Bins with a 
number of counts N - have been set to log = -1 and a constant equal to 7.0, 4.5 and 2.0 has been added to logA^ for ad, be 
and cf, respectively. 
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Fig. 8. a Colour-magnitude diagram showing the ZAMS locus 
(thin solid line), three isochrones (10^, 10'' and 10^ yr) from 
the pms tracks by Palla & Stahler (1999; dashed lines) and a 
"mean" locus used to deredden NIR sources (thick solid line), 
all for a distance of 700 pc. A few spectral types on the ZAMS 
and star masses on the 10^ yr isochrone are labelled, along with 
the ages of isochrones. An arrow shows the reddening for A v = 
10 mag according to Rieke & Lebofsky (1985). The error made 
in dereddening a O.l-M© 10^-yrs-old pms star is indicated by 
the horizontal short-dashed lines, b same as a, but only the 
ZAMS locus, the 10^ yr isochrones from Palla & Stahler (1999; 
labelled as PS98) and D'Antona & Mazzitelli (1994; labelled 
as DM94) and the mean locus are drawn, c same as a, but only 
the ZAMS locus, the lO*" yr isochrone from Palla & Stahler (at 
a distance of 500, 700 and 900 pc, labelled) and the mean locus 
are drawn. 

ing from 10^^ to 10** years in steps of 0.5 (in logarithm). From 
these, we got the actual "unreddened" KLFs with and without 
a random NIR excess component. In both cases, a variable red- 
dening was added to each synthetic member and the KLFs were 
then dereddened, using both the original mean lous and the one 
shifted in H - K hy 0.4 mag to compensate for the NIR ex- 
cess. Each of the "dereddened" KLFs was compared to all the 
"unreddened" ones and the mean squared differences (x^) in 
counts per magnitude bin (of width 0.5 mag) were computed. 
The results can be summarized as follows: 

1 . when no NIR excess is included, the least is obtained for 
the same age; 

2. when a NIR excess component is included, the best fit is 
generally obtained for the same age, but an ambiguity arises 



between the 10^'^ yrs old and the 10^ yrs old models. The 
same occurs when shifting the mean locus m H - K, al- 
though the obtained values are lower than before; 

3. using the evolutionary tracks of D'Antona & Mazzitelli 
(1994) in synthesizing the clusters, the best fits are always 
obtained for the same age, excepted that the dereddened 
KLF for a 10^ yrs old cluster (including a NIR excess com- 
ponent) is equally well fitted by unreddened KLFs 10^'^ and 
10^ yrs old; 

4. when the cluster distance is set to 500 pc, the best fits are 
for the same age, but when NIR excess is included, the 10^ 
and 10^ yrs old "dereddened" KLFs look younger by 0.5 (in 
logarithm). Shifting the mean locus m H - K corrects for 
the effect, but causes the 10^'^ yrs old "dereddened" KLF 
to be fitted by the 10** yrs old KLF; 

5. when the cluster distance is set to 900 pc, the best fits are 
always for the same age, excepted that shifting the mean 
locus mH - K causes the 10^^ yrs old "dereddened" KLF 
to be well fitted by the 10^ yrs old KLF. 

We performed similar tests by constructing 3 new syn- 
thetic clusters with continuous star formation in the ranges 
^05.5 _ iq6.s^ _ YQi.i and 10^-'* - lO'*-^ yrs. Now, the "dered- 
dened" KLF is always best fitted by the "unreddened" KLF for 
the same cluster age range, unrespective of the presence of a 
random NIR excess. Comparing the "dereddened" coeval KLFs 
to the 3 unreddened ones with continuous star formation, indi- 
cates that only for the oldest coeval clusters some ambiguity 
may arise with KLFs of evolved cluster with continuous star 
formation. In conclusion, our dereddening technique allows to 
reconstruct the actual KLF with sufficient accuracy provided 
that the adopted mass-luminosity relation and distance to the 
VMR are within reasonable values and that no NIR excess is in- 
cluded. When the latter is not negligible, the dereddening tech- 
nique stretches the "dereddened" KLF at the high luminosity 
end with respect to the "unreddened" one, but the age of the 
cluster appears still discernible within a factor of 10*^^. 

By dereddening the KLFs of the observed clusters as ex- 
plained, we found only small differences between the dered- 
dened KLFs either using the mean locus or the shifted one. So, 
at last we decided to deredden the clusters' KLFs onto the un- 
shifted mean locus. This also because, as can be seen in Fig.|3 
only a fraction of sources exhibit a NIR excess, whereas our 
tests always assumed that all sources had a NIR excess compo- 
nent. Furthermore, it is likely that the most prominent NIR ex- 
cess is displayed by the brightest sources, i. e., those falling 
in a luminosity interval where the KLFs are anyway scarcely 
sampled due to the poor statistic. Thus, the dereddened KLFs 
should be well representative of the actual unreddened ones. 
These are shown in Fig. It is evident that all are quite 
similar, both in extent and in shape. An estimate of the com- 
pleteness limit has been computed as follows. A vertical line is 
drawn enclosing the maximum H - K displayed by the sources 
in the CMDs and its intersection is found with the line of the 
photometric completeness limit. Projecting this point along the 
reddening vector onto the mean locus then yields the "dered- 
dened" completeness limit. This is estimated to be ~ 14 and 
it is shown in Fig.[nas a vertical dotted line. 
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Table 3. Significance level of the hypothesis that two cluster 
KLFs are the same distribution from the chi-square test. The 
crossing of each line and column lists the significance level for 
that pair of KLFs. 



Field 
name 


IRS 16 


IRS 17 


IRS 18 


IRS 19 


IRS20 


IRS21 


IRS 16 




0.95 


0.93 


0.26 


0.32 


0.79 


IRS 17 


0.95 




0.61 


0.06 


0.52 


0.25 


IRS 18 


0.93 


0.61 




0.20 


0.34 


0.71 


IRS 19 


0.26 


0.06 


0.20 




0.04 


0.96 


IRS20 


0.32 


0.52 


0.34 


0.04 




0.28 


IRS21 


0.79 


0.25 


0.71 


0.96 


0.28 





All KLFs appear extremely flat across the "dereddened" 
completeness limit, which corresponds to pms stars of ~ 0.1 
Mo for an age of 10^ yrs (~ 0.4 Mq for an age of 6 x 10* 
yrs), according to the tracks of Palla & Stahler (1999). This 
indicates that the underlying IMF is flat or slightly increas- 
ing across this mass. A break on the slope of the KLFs can be 
roughly identified at > 10, suggesting a break in the IMFs 
at < 2.5 M0 for 10* yrs (< 3.5 Mq for 6 x 10* yrs). Above this 
mass limit, the source statistic is quite poor The magnitude 
shifts of the KLFs with respect to each other due to small dif- 
ferences in the clusters' distance can be estimated as follows. 
The maximum projected distance between clusters is ~ 1.5°, 
corresponding to ~ 18 pc. Hence, the maximum distance dif- 
ference along the line of sight is expected of the same order. 
This would cause maximum variations in the distance modulus 
~ 0.06 mag, much less than the binning width. 

The hypothesis that all KLFs arise from a same distribution 
can be tested now. If plausible, this would confirm that both the 
underlying IMF and the star formation history can have been 
similar for all clusters. Hence, we performed the chi-square test 
on all pairs of KLFs; the obtained significance levels are listed 
in Tab. 13 The KLFs are compared in the range Ks - 10 - 14; 
the lower limit is chosen to avoid counting problems due to 
saturation of the array. It can be seen that the hypothesis can- 
not be disproved at generally high significance levels. IRS 16, 
17 and 18 exhibit quite high probabilities to represent a same 
distribution, whereas IRS 19 and 20 are less likely to meet the 
hypothesis, but still this cannot be disproved. Performing the 
test on the dereddened KLFs uncorrected for the field contribu- 
tion generally causes slightly higher levels for the larger signifi- 
cance cases and lower levels for the smaller significance cases, 
but of the same order as before. In summary, the probability 
that the observed KLFs are drawn from a same distribution ap- 
pears significant. 

4. 7. Deriving tine initial mass function of the clusters 

As a last step, we want to derive the IMF underlying the clus- 
ters' KLFs. So far, most of the efforts in the literature have been 
devoted towards constructing synthetic KLFs from a given IMF 
and star formation history (see, e. g., Lada & Lada 2003). These 
are then fitted to the observed KLFs and the IMF underlying 
the best-fit KLF is thus extracted. Due to the hmited statistic of 



I 



10 - 




6 8 10 12 14 16 18 

K (mag) 

Fig. 11. The KLFs of the 6 clusters, corrected for the field stars 
contribution and dereddened according to the method described 
in the text. The vertical dotted line marks the estimated com- 
pleteness limit. 

our clusters, we chose to directly obtain the IMF by using the 
relation: 

dN ^ djV dlogM 1 
dlogM dK^ dK ' ' ^ ^ 

Note that this IMF is still model-dependent since the mass- 
magnitude relation has to be derived from a set of evolutionary 
tracks. And a star formation history has to be assumed. 

We computed the mass-magnitude relation from the evolu- 
tionary tracks of Palla & Stahler (1999) for a grid of ages span- 
ning 10^ to 10^ yrs. If the clusters were coeval, Eq.[0 could 
be applied straightforwardly and the IMF could be obtained di- 
rectly from the KLF. This because the freezing of time allows 
one to fix the mass-magnitude relation and a given magnitude 
corresponds uniquely to a given mass (excepted in small mag- 
nitude intervals; see below). However, in order to apply Eq.^ 
to the clusters' KLFs, both the KLFs and the derivative of the 
mass-magnitude relation have also to change slowly. Whereas 
this is true for the KLFs, it can be seen in Fig.^jthat there are 
limited intervals in magnitude where this condition is not sat- 
isfied by the derivative of the mass-magnitude relation. In this 
respect, there appear to be two kinds of problems. The first one 
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Fig. 12. Mass-magnitude relation for pms stars 10^ (plusses), 
10* (asterisks) and 10^ (open diamonds) yrs old, from the evo- 
lutionary tracks of Palla & Stahler (1999). Masses range from 
6 (upper end) to 0. 1 (lower end) M0 (a few points are indi- 
cated). Fits to the tracks are shown as dashed lines. The mass- 
magnitude relation for ZAMS stars is also drawn (M > 6 Mq.). 



is evident for very young ages (see the track 10'' yrs old), where 
dlogM/d^T — > 00 around 1.5Mq. This is a signature of the 
birthline (see, e. g.. Fig 10 of Palla & Stahler 1993) and, hence, 
reflects the initial conditions of the pms contraction. The other 
problem is common to all ages and is represented by a sort of 
"tooth" moving from higher to lower masses at increasing ages. 
This causes a small magnitude interval to be "degenerate" in 
mass (i. e., stars with different masses display the same K). This 
appears to be related to the end of the radiative contraction and 
the start of the main sequence. We avoided an interval spanning 
1 magnitude around the points where dlogM/d/T — > 00 and 
the whole range (increased by 0.25 mag at each end) where the 
masses are "degenerate". Outside these intervals, we used poly- 
nomial fits to approximate the mass-magnitude relation (see 
Fig. \121 . In turn, these have been used in order to convert K 
into masses and to estimate d log M/dK. 

As for massive stars, according to the Palla & Stahler sce- 
nario an object enters the ZAMS after accreting 8-10 Mq. 
Hence, we used the mass-luminosity relation as derived from 
the ZAMS for stars with M > 6 Mq. We adopted the abso- 
lute magnitudes given by Panagia (1973) and the colours by 
Koornneef (1983) to obtain the K magnitudes of OB stars. The 
masses were then determined by converting the bolometric lu- 
minosities given by Panagia (1973) according to the formulae 
by Tout et al. (1996). This is also shown in Fig.^] For an as- 
sessment of the errors affecting this relation see, e. g., Vacca et 
al. (1996). 



The clusters' IMFs for assumed ages of 10*, 3 x 10'' and 6 x 
10* yrs and coeval star formation, are shown in Fig.[Ol These 
are compared with the IMFs from Scalo (1998) and Kroupa et 
al. (1993). The Scalo's IMF (in the form dN/dlogM) is ap- 
proximated by a power law with F = -1.7 between 10-1 M0 
and F = -0.2 between 1-0.1 M0. It is derived by a compila- 
tion of different determinations. That of Kroupa et al. (1993) is 
representative of field stars and is similar to the Scalo's (1998). 
In particular, F = - 1 .7 for M > 1 Mq, F = - 1 .2 between 1-0.5 
Mq and F = -0.3 between 0.5-0.08 Mq. 

The errors, propagated from the statistical uncertainties of 
the KLFs, are drawn only for 10* yrs old. Because of the in- 
tervals to be avoided, the IMF 10* yrs old is not defined for 
M > 1.6 Mq whereas the one 3 x 10* yrs old exhibits a gap 
between 1.6 and 4 Mq. Moreover, at this stage we do not show 
the IMF for the most massive stars through the ZAMS mass- 
luminosity relation because of the poor statistic. However, it 
is clear that the clusters' IMFs are consistent with the "com- 
parison" ones within the uncertainties. In the case of IRS 19, 
the shape of the IMF seems to be aff'ected by the poor statistic 
rather than by actual differences with the comparison IMFs. 

In a few cases, a better agreement would be obtained by 
shifting the breakpoint envisaged in the Scalo (1998) IMF from 
1 to ~ 1 .5 Mq, especially if older ages are assumed. But it is 
difficult to assess whether this reflects an actual difference. We 
also checked that using the evolutionary tracks of D'Antona 
& Mazzitelli (1994) does not affect significantly the obtained 
IMFs, excepted for 10^ yrs. 

Assuming a coeval star formation is obviously an oversim- 
plification. Nevertheless, if star formation is peaked at a given 
age, this should not change the results. Furthermore, it is ev- 
ident from Fig. El that the derivative of the mass-magnitude 
relation changes slowly with time between 1-2 and 0.1 Mq. 
What changes is the K magnitude. We can hope that the re- 
sults should be still little affected even in the case of continu- 
ous star formation if it lasts for a time such that the change in 
K is ^ 0.5 mag (i. e., the binning interval of our KLFs). This 
probably translates, e. g., into a duration of ^ 2 x 10* yrs for 
ages of a few 10* yrs. Some concerns arose from the possibil- 
ity that star formation may span ~ 10^ yrs and accelerate, as 
proposed by Palla & Stahler (1999). Hence, we approximated 
clusters of ~ 100 - 200 members with accelerated star forma- 
tion by adding a KLF synthesized in the coeval 10* yrs old case 
to a KLF synthesized in the case of continuum star formation 
from 10^ to 10* yrs, with different fractions of the total mem- 
ber number. Then, we derived the IMF from the cumulative 
KLF by using our method. We checked that the obtained IMF 
is very similar to the underlying one (Scalo-like) in the 1.5-0.1 
Mq range even if up to a 35 % of the synthesized population is 
produced through continuum star formation. In particular, in a 
log(dA^/d log M) vs. logM plot the purely coeval IMF and the 
mixed one still coincide within the statistical errors. This may 
be understood analytically: the continuum component of the 
KLF is flatter and spread on a larger magnitude interval than 
the coeval component, and its underlying area is also much 
smaller. But we checked that even increasing the continuum 
fraction of the synthesized population to a 50 %, the obtained 
IMF does not change dramatically, although appearing slightly 
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flatter. However, this assumes that the underlying IMF is the 
same in the two population components; only in the unlikely 
case that stars of different masses are originated at very differ- 
ent times we could expect major differences between the ob- 
tained and the actual IMF in the 1.5-0.1 Mq range. 

We could reverse our point by trying to constrain the clus- 
ters' age assuming that the underlying IMF is a standard one, 
like our comparison IMFs. If so, the clusters' KLFs are consis- 
tent with coeval star formation and an age between 1-6x10^ 
yrs. Younger ages can be discarded since they produce IMFs 
peaked at 0.3-0.1 Mq with Salpeter-like slopes. Also, clusters 
older than 10^ yrs can be discarded, since the resulting IMFs 
would be peaked at > 2 Mq and declining towards higher and 
smaller masses. 

5. Discussion 

Knowing the dereddened cluster KLFs allows us to estimate 
the number of members A^r up to the completeness magnitude 
(i. e., down to 0.1-0.4 Mq), to be compared with Ic in Table|2l 
The results are listed in Table 0] and it can be noted that A^r is 
roughly twice the corresponding value of Ic. Clearly, the latter 
underestimates the true number of members, probably due to 
the fact that the sky contribution is estimated in a region still 
contaminated by stars from the cluster population. In fact, de- 
riving the clusters' KLFs within a radius of 1 arcmin (instead 
of 2 arcmin) results in a much lower Nj. 

If we want to include in A^r all clusters' members down to 
the hydrogen burning limit, then the values derived from our 
observations have to be corrected for the age. Having proved 
that the clusters' IMF is very similar to standard ones, like those 
of Scalo (1998) or Kroupa et al. (1993), by adopting for it the 
IMF of Scalo (1998) we can estimate the correction to apply 
to A^r. E. g., for a 6 x 10^ yrs old cluster, our completeness 
magnitude corresponds to a mass limit of 0.4 Mq. Integrating 
the Scalo IMF from 0.4 to 0.1 Mq yields ~ 55 % of stars in 
that mass range, so the observed A^j should be increased by 
a factor of 2.2. This largely exceeds the statistical (counting) 
uncertainties listed in Table|2] Note that a further underestimate 
is likely to arise from having selected sky areas close to the 
clusters in order to construct the control KLFs (i. e., possibly 
still containing a small number of clusters' members). And we 
do not consider here the brown dwarf population (which is not 
sampled by our observation). 

Actually, A^r is bound to be underestimated anyway, due to 
the unknown fraction of close companions which escape obser- 
vations. How critical this problem may be is demonstrated by 
Fig. 6 of Giannini et al. (2005); using ISAAC images of IRS 
17, these authors can resolve source # 40 of Massi et al. (1999) 
into 6 different components. We checked that the same star in 
our Ks image is resolved into only 3 components, although the 
brightest ones. 

The total mass (down to the hydrogen burning limit) of the 
clusters can be derived, as well, by integrating the Scalo IMF 
times the mass. The results are listed in Table |3 and, like A^x, 
can be underestimated by a factor 2.2 or so, depending on the 
age. The highest expected stellar mass x in each cluster was 
estimated by integrating the Scalo IMF, as well, from oo to a 



mass X and requiring that the result is 1 . The obtained values 
(see Table ^ are consistent with the bolometric luminosities 
given by Massi et al. (2000, 2003) and confirm that the most 
massive stars in the 6 clusters are intermediate-mass stars. Note 
that an increase in A^r of a factor 2.2 would roughly double x. 

At last, we can derive the star formation efficiencies of the 
parental clumps if the mass of the associated gas is known. 
We used the mm (C'^O) observations of Wouterloot & Brand 
(1999) and Yamaguchi et al. (1999). Because of the low resolu- 
tion of the latter data, these probably refer to matter on a larger 
scale than found by Wouterloot & Brand (1999), i. e. clouds 
rather than clumps or cores, and are likely to underestimate 
the gas content of these larger scale structures. Nevertheless, 
all derived efficiencies are < 26 %, excepted towards IRS 20 
where a star formation efficiency of ~ 80 % is found. These 
are typical values for young embedded clusters (Lada & Lada 
2003). Instead, if taken at face value, the high star formation 
efficiency towards IRS 20 implies that the cluster is emerging 
as a bound structure (Lada & Lada 2003); however, we remark 
that the area around IRS 20 covered by the mm observations 
of Wouterloot & Brand (1999) is very limited and the gas mass 
could then have been underestimated. 

Although our NIR observations were able to clarify the 
shape of the clusters' IMF between 1-3 and 0.1-0.4 Mq, they 
fail to strongly constrain the clusters' age. This would have al- 
lowed to more accurately model the dereddened KLFs and de- 
rive the underlying IMF. Also, the substellar mass regime is not 
probed by our observations, requiring much deeper imaging. 
However, as shown in Fig. ^2 the dereddened KLFs maintain 
roughly the same slope above the completeness magnitude (re- 
member that 5 fields are ~ 1 mag deeper than the 6th one). 
No strong decrease is evidenced; hence, because of uncom- 
pleteness, the KLFs can at most increase above this limit. This 
means either that the clusters are slightly older than 1 x 10^ 
yrs and stars fainter than Ks = 14 are still above the hydrogen 
burning limit (and follow a Scalo-like IMF), or that the clus- 
ters are as old as, or younger than, 1 x 10*' yrs and the IMF 
itself is still increasing at the hydrogen burning limit (assumed 
dlogM/d/T does not exhibit major changes below 0.1 Mq). 

As for the high mass end of the IMF, this is not only affected 
by the poor statistic. The brightest sources in our JHKs images 
are also saturated and, probably, most of the brighter sources 
are largely affected by a NIR excess emission. However, the 
data obtained by Massi et al. (1999, 2003) are less affected by 
saturation or non-linearity problems. As for IRS 16, it coin- 
cides with an HII region and these authors found that it is ex- 
cited by a B0-B5 source, which is the most massive cluster's 
member (consistently with the highest expected stellar mass 
listed in Table|4] given also the uncertainty on A^r)- In the case 
of the remaining clusters, the most massive clusters' members 
are the sources emerging in the NIR and identified by Massi 
et al. (1999) as the counterparts of the IRAS point sources. 
Whereas the clusters' IMF is consistent with the mass of the 
most massive star in each cluster, imaging surveys of IRAS 
point sources in the VMR have so far failed to unveil high 
luminosity IRAS sources associated with isolated young stars 
(Massi et al. 1999, 2003). 
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Table 4. Clusters' properties derived from their IMF. Column 
2 lists the number of cluster members (down to the mass com- 
pleteness limit), Col. 3 the cluster total mass (as in Col. 2), Col. 

4 the mass of the most massive star within the cluster and Col. 

5 the star formation efficiency. 



Field 


Members 


Mass 


Largest 


SFE 


name 






star 








(Mo) 


(Mq) 




IRS 16 


144 ± 15 


103 


7.1 


0.24" 


IRS 17 


126 ± 15 


90 


6.5 


0.26* 


IRS 18 


124 ± 14 


89 


6.5 




IRS 19 


69± 13 


49.5 


4.4 


0.095* 










0.21° 


IRS 20 


98 + 14 


70 


5.5 


0.8* 


IRS 21 


95 ± 14 


68 


5.4 


0.11° 



° gas mass from Yamaguchi et al. (1999) 
* gas mass from Wouterloot & Brand (1999) 



A way to investigate the high mass end of the clusters' IMF 
is by adding together all 6 IMFs, as to increase the statistic. 
This should also approximate the field star IMF, if most stars 
originate in clusters and then disperse. The "global" IMF is 
shown in Fig. The four most massive bins have been ob- 
tained by using the ZAMS mass-magnitude relation discussed 
in Sect. 14.71 We dropped three points with resulting masses 
> 20 M0. We checked that these are due to the young stars 
or protostars lying in the centre of IRS 17, IRS 19 and IRS 
20 and identified by Massi et al. (1999), then members of the 
clusters. But, as shown by Massi et al. (2000, 2003), based 
on the measured bolometric luminosity and the upper limit on 
the radio continuum emission, they are intermediate-mass, not 
high-mass (proto-) stars. Clearly, they exhibit a strong NIR ex- 
cess emission and the ZAMS mass-magnitude relation cannot 
be applied to them. 

The global IMF is quite similar to the comparison ones, 
particularly at > 10^ yrs. The IMF 10'' yrs old implies a break 
at M > 1 Mq, but the mass of many objects in the intermediate- 
mass regime may have been overestimated because of the NIR 
excess. At the high-mass end, the IMF slope appears slightly 
steeper than that of the comparison ones, but here the statistic 
is still poor However, it is noteworthy a drop-off of the IMF 
at M ~ 10 Mq. There are 13 stars in the 4 most massive bins 
(M > 6 Mq), of which 1 in the most massive one (M ~ 14 
Mq) and 2 in the next most massive one (M ~ 10 Mq). This 
indicates that none of the 6 clusters could produce a high-mass 
stars, although they all originated intermediate-mass stars as 
predicted, e. g., by the Scalo's IMF based on their number of 
members. But applying the same statistical argument, since we 
are now sampling a total number of stars > 656 (see Table|^ we 
should find a star with M > 22.5 Mq (i. e., an 08-09 ZAMS 
star). Such an object would have a bolometric luminosity > 10^ 
Lq and form an HII region. There is not such a massive source 
within the 6 clusters (Massi et al. 2000, 2003). Hence, the drop- 
off appears to be a real feature of the global IMF. 

If the actual clusters' IMF is a standard one, then the most 
massive stars in each cluster depend on the cluster mass, as 
expected based on statistical grounds. But they fail to follow 
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Fig. 14. Sum of the 6 initial mass functions, same symbols as 
Fig. The four most massive bins have been derived using 
the ZAMS mass-magnitude relation discussed in the text. 

statistics when added together, suggesting that the dependence 
on the cluster mass is "physical". I. e., high-mass stars need 
large clusters to be formed. On the other hand, it may be that 
the IMF in the whole region is a standard one, while that of the 
single clusters is much steeper at the high mass end. Then, we 
should find massive stars elsewhere in the region. Massi et al. 
(1999, 2003) do not find massive stars in their D-cloud sample 
of IRAS sources, but those objects might well have escaped the 
selection criteria adopted by Liseau et al. (1991). This would 
reconcile with a statistical view of the IMF, but if the clusters' 
IMF is steeper than that of field stars or open clusters, it is 
likely so because of different physical conditions in the parental 
gas, so, in a sense, still "physical". A careful multi-wavelength 
examination of the whole region is then necessary in order to 
settle this point. 

6. Conclusions 

We have presented the results of a Near-Infrared photometric 
survey of a sample of 6 embedded young clusters, all associated 
with luminous IRAS sources in the Vela D molecular cloud. 
They lie at the same distance and have been imaged with the 
same instrumental setting; this allowed us to derive the clus- 
ters' properties in a consistent way, without most of the obser- 
vational biases often affecting the comparison of different star 
forming regions. The clusters' K Luminosity Functions have 
been obtained from the same images, by subtracting a control 
K Luminosity Function, for sources > 2' from the clusters' 
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centre, from that for sources < 2' from the clusters' centre. A 
simple method for dereddening the K Luminosity Functions, 
based on colour-magnitude diagrams, has been developed and 
described. The main results are the following: 

1. Our images confirm the previous results of Massi et al. 
(2000, 2003), i. e., all selected fields contain young em- 
bedded clusters with a radius ~ 0.1 - 0.2 pc and a mean 
stellar density ~ 10^ stars pc"^. We also resolved, for the 
first time, IRS 17 into 2 sub-clusters. 

2. the dereddened clusters' K Luminosity Functions have 
been compared in pairs using the chi-square test; the hy- 
pothesis that all arise from a same distribution cannot be 
disproved at significance levels ranging from 0.04 to 0.95 
(but > 0.50 in most cases). This suggests that the under- 
lying Initial Mass Function and star formation history may 
have been quite similar. 

3. The clusters' Initial Mass Functions, obtained from the 
dereddened K Luminosity Functions in the hypothesis of 
coeval star formation, are compatible with standard ones, 
like those of Scalo (1998) and Kroupa et al. (1993), for 
assumed ages ranging from 10^ to 6 x 10^ yrs. Major de- 
partures would occur if the clusters were either younger or 
older than assumed. However, they are unUkely to be older 
than 5 X 10* yrs since all of them are still associated with 
molecular gas. 

4. The global IMF, obtained by adding together all 6 single 
IMFs, is similar to a standard IMF. But it appears steeper 
at the high-mass end and exhibits a fall off" at ~ 10 M0. 
E. g., on statistical grounds the Scalo IMF would predict a 
star of > 22.5 Mq within one of the cluster, which is not 
found. Then, either high-mass stars need larger clusters to 
be formed, or the single IMFs are steeper at the high-mass 
end because of the physical conditions in the parental gas. 

5. Based on their K Luminosity Functions, the clusters' mem- 
bers range from 69 (IRS 19) to 144 (IRS 16). Assuming 
an Initial Mass Function like Scalo (1998) yields total clus- 
ters' masses ranging from 50 to 100 M©. If the clusters are 
10* yrs old, the stellar mass is sampled down to 0.1 M©. If 
the clusters are as old as 6 x 10* yrs old, both number of 
members and mass have to be increased by a factor ~ 2.2 
to reach the same mass lower limit. 

6. The star formation efficiency ranges from 0.1 to 0.8, being 
< 0.26 in all cases but one (IRS 20). 
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Fig. 13. The initial mass functions for the 6 clusters, derived assuming coeval star formation and an age of 10^ (open squares), 
3x10^ (asterisks) and 6x10^ (open diamonds) yrs. All cluster IMFs are compared with the IMFs from Scalo (1998; dotted Une) 
and Kroupa et al. (1993; dashed Unes), each with two different normaUzations. For the sake of clarity, errorbars are drawn only 
for lO*" yrs old. 



